arXiv: 1508.06146v2 [cond-mat.soft] 24 Nov 2015 


Wrinkling Reveals a New Isometry of 
Pressurized Elastic Shells: Supplementary 

Information 

By Dominic Vella, Hamid Ebrahimi, Ashkan Vaziri and Benny 

Davidovitch 


1. Introduction 


In this Supplementary Information, we present more details of the calculation pre¬ 
sented in the main text. In particular in ^we present the details of the calculation 
of the base shape, Co(r), about which wrinkles occur. The main results of this anal¬ 
ysis (which are also given in the paper) are summarized in In (21 we estimate 
the energy of the wrinkle pattern and thereby determine the scaling law for the 
wrinkle number given in the main paper. Finally, in 21 we compare the results for 
the indentation of a pre ssurized she l l with those obtained for the indentation of an 
indented floating sheet ( Vella et al . 2015ll . 


2. Governing equations 

We begin with the dimensionless equations given bv IVella et al. ( 2011 ) for the 
deformation of a pressurized shell. For consistency with the main text we use (^(r) 
to denote the normal displacement of the shell from its pressurized, pre-indented 
state. We set 
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( 2 . 1 ) 
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where Y = Et is the stretching stiffness of the shell. With this non-dimensionalizatio n, 


we fin d that the shallow shell equations (|Ventsel fc Krauthammeii . l200ll : IVella et al. 
20111) become 
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(rH)i/2 

is a dimensionless measure of the tension within the shell due to inflation alone, 
CToo = pR/2. 
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In the main text, we focus on the membrane limit, which corresponds to r = 
oo. We shall therefore neglect the term proportional to in (12.21) . However, 

when discussing the scaling behaviour of the number of wrinkles, we shall return 
to consider the role of this term. 

A key control parameter is the indentation depth, which in these dimensionless 
variables is 

5 = SR/ll- (2.5) 

This imposed indentation depth is imposed by the boundary condition Z{0) = —<5, 
which in turn requires a point force F applied at p = 0. A key quantity of interest 
is the relationship between F and <5. 


(a) Tensile regions 

In regions where both the radial and hoop stresses remain tensile, i.e. tr^rj crge > 
0, the shell deformation is described by the membrane-shell equations 


and 


lA 
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ipd>) 


p dp V dp y 
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^ dp 2 V dp y 


( 2 . 6 ) 


(2.7) 


(6) The wrinkled region 

It has been shown previously IVella et al\ ( 2011 1 that the hoop stress becomes 
compressive {age < 0) for sufficiently large indentation depths, 5 > 5c ~ 2.52. In 
reality, a membrane cannot support a compressive stress and so we must instead 
assume that the compression is completely relaxed, agg 0, within an, as yet 
unknown, wrinkled region Lj < p < Lq- In this region, the equilibrium equation 


dfj 

rr ^ (Jrr ^90 

dr r 


( 2 . 8 ) 


leads immediately to 

C 

arr = —■ (2-9) 

r 

Equation (12.91) replaces the equation of compatibility of strains (12.71) in wrinkled 
regions. 

The appropriate expression of vertical force balance in a wrinkled region may 
be found by substituting a dimensionless version of (12.91) . oi = C' = C|pR^p, 
into (EH)- Dropping the primes, we have 


dp2 “ 
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(c) Matching boundary conditions 

Equations were solved in the tensile regions Q < p < Li and p > Lq 

while equation (12.101) was solved for the wrinkled region Li < p < Lq- The complete 


set of boundary conditions applied were 
Z(0) = —(5, lim \p^' — 

[Z]% = [Zf^^ = [Ur]% = 0 , 

[zf^ = [Zf^ = [Urf.^ = 0, 

Lq 

Z —>■ 0, 'k ~ p/2, as , 


= 0 


(2.11) 

vE-(L7)=0, 

VI/'(L7) = 0 

(2.12) 

A 

A 

II 

Q 

vI-'(L+) = 0 

(2.13) 

) —>■ oo. 


(2.14) 


3. Analytical results in the limit of large indentations 

From the numerics, it appears that L/ —>■ 0 as 5 increases. At the same time, we 
expect the general horizontal scale to be increasing (specifically ^ 5'^/'^) and so we 
seek to exploit the fact that Li/d^l'^ <C 1 to simplify the problem within the inner 
tensile region, p < Li', in particular, we will be able to neglect the effect of the 
pressure loading within this region. 


(a) The inner tensile region 0 < p < Lj 

The numerical solution of the reduced, FT, problem suggest that Lj ^ 5“^/^ 
and C ^ (5^/^ for 5 ^ 1. We therefore let L/ = , C = and rescale the 

equations by letting 


e = p/L/, io{f) = Z/l X = VI//C = vI//^^i/2. 

We find that 

i A _ Ai A = A _ 

f dC f df de J p~s 2 ttXip f ■ 


(3.1) 


(3.2) 


While it may seem that there is nothing that can balance the second term on the 
LHS of (13.21) . it is important to remember that we do not expect dw/d^ to be 
0(1) within the tensile region; rather this quantity will be small, and so the correct 
balance is between the second term on the LHS and the second term on the RHS. 
Integrating, we find that 


dw FXi 
^ d^ 27r/i5^ 


= F. 


(3.3) 


This now gives us an expression with which to eliminate uj' from the equation of 
compatibility, which in rescaled coordinates reads 


4 




Xi duj 16^ 


(3.4) 


Assuming that uj' > 0{S ^), i.e. that F > 0(1), the first term on the RHS, which 
arises from the curvature of the shell, may be neglected in comparison to the second, 
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which arises from the strain caused by out-of-plane deformation. We then find that 


4 


4(?x) 


L? dC 


15^ / dw 


2 A//i \ d^ 


1(52J-2 


-2 


2 A//i 


(3.5) 


Eqn (13.511 is the usual compatibility of strains equation that is found in the FvK 
equation for a naturally flat sheet: the curvature of the shell does not play a role 
in the inner tensile region. 

Equation (13.511 is to be solved on 0 < ^ < 1 with boundary conditions 


arr{Li)=ClLi 

(yee{Li) = 0 


x(i) = 1 
x'(i) = 0 


(3.6) 

(3.7) 


and zero horizontal displacement at the origin, i.e. limp_^o \^P^' ~ = 0 . 

Eqn (13.51) with boundary conditions (I3.6I) - (I3.7I) can be so lved analytica l ly us¬ 
ing steps s i milar to those outlined in related problems ('see IChooin et al,. . 2008 : 
Vella et aJ.. . 20151 for example). We first let: 


f = 


so that 


A/^ ’ 

d 2 $ 

dr]'^ 


V = ^ , $ = 


Jl. 

8$2 


with boundary conditions 

$( 0 ) = 0 , $'(!) = 1 / 2 , $( 1 ) = 1 . 

The first integral gives 


d$ 

d?7 


/ fl + Ad>Y^^ 
2 


d> 


- 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


where (1 -I- = / ^ using the boundary conditions at 77 = 1. Integrating again, 

we have 


A 


(1-r/) = sinh-i(yl$)i /2 _ ^ 1/2 

(3.12) 


2(1-P .4)1/2 

The condition at 77 = 0 then gives that the constant A must satisfy 

= (l+yl)i/2_yi-i/2sinh-iyli/2. 


A 


2(1 +Ay/2 


(3.13) 


which has solution A fa —0.696529 and hence / ~ 1.81527. Note that this is pre¬ 
cisely the sa me result as was fo und for the indentation of a floating, naturally flat, 
membrane bv I Vella et al. ( 2015 1. This similarity occurs because sufficiently close to 
the origin neither the internal pressure nor the curvature of the shell play a role: 
the balance is purely that between the indentation force F and the tension within 
the deformed membrane. 
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The result / « 1.81527 shows that IF ~ 5“^ and hence the dimensionless inden¬ 
tation force F ^ FS'^ ^ 6. 

However, to make further progress we need to understand the behaviour of the 
solution in this region close to the outer region of its validity (i.e. close to the inner 
radius of the wrinkle pattern, p = L/, = 1). We note that (13.31) can be rearranged 

to give 

from which we find 

o j: 

^ = -l + ]^sinh-i(H$)i/2. (3.15) 


Crucially, we see that the rescaled vertical deflection at the inner portion of the 
wrinkled region = 1) = —1 -|- 0{l/5)\ the boundary condition for the wrinkled 
problem at p = Lj should, at leading order, be simply Z = —6 — the imposed 
indentation depth. This is the result used in the simplified calculation presented in 
the main text. 

As well as matching vertical displacements between regions, we must also match 
the slopes of the displacements (as a simple force balance argument illustrates). We 
therefore also need to know the value of at p = LJ , which we find from ()3.3I) to 
be 


dp 


= Fp/^/Xi = —^ 

p=L7 


2-kC 


(3.16) 


(6) The wrinkled region: Lj < p < Lq 

In the wrinkled region we have 'k = C and that the vertical displacement satisfies 
(I2.10L which may be integrated immediately to give 

Z = p^/2-^+A + Bp. (3.17) 

The constants A and B can be determined from the continuity conditions at the 
inner boundary, p = Lj, where we have 


Z{Lj) 


-I 


dZ 


dp 


p=Li 


F 


(3.18) 


Substituting these conditions into (13.171) and noting that L/ —>■ 0 as 5 —?► oo imme¬ 
diately gives A = —S, B = F/(2TrC) (the latter finding being consistent with the 
appearance of a Dirac ^-function in (12.101) 1. Hence, we have 


We shall proceed to patch this profile onto that outside the wrinkled region, 
p > Lo- In this regard, it is useful to note that differentiating the profile (13.191) and 
evaluating at p = Lq gives 


dZ 


dp 


p=Lo 


F 


+ Lo 



(3.20) 
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(c) Beyond the wrinkled region: p > Lq 


Eqn (12.61) can be integrated once to give 

dZ A + p^/2 

The constant of integration A is immediately found by continuity of Z' and 'I' at 
p = Lq to be ^ = —FI^-k. (Again, this choice of the constant of integration is 
reassuring since it implies that the force being supported in the far field is that due 
to indentation, F.) We therefore have 


dZ F/tt — p^ 

d^ 2^ 


(3.22) 


To solve the problem completely, however, we must determine the Airy stress 
function 4'. To do this, we eliminate Z' from the equation of compatibility (EJl) 
using (13.221) to give 


dp 


1 d 

pdp 


(P^) 





(3.23) 


The solution of eqn (13.231) must match continuously onto that in the wrinkled region 
and so we must have that 


M/(L+) = C, vI;'(L+) = 0. 


(3.24) 


(Note that far away from the wrinkled region, we expect that ^ p/2, which will 
emerge necessarily from the structure of the equation.) 

To make progress, we again make the substitution used in the inner tensile 
region, i.e. $ = p'l', p = p^, with which (13.231) becomes 


d^<i) 1 / F/tt — ryX 

d^ ~ 8 V 2$ J ’ 


(3.25) 


with boundary conditions $(Aq(5) = 4>'(Aq 5) = p/(2Ao) where Xq = 

The condition at oo becomes that $ ~ 77/2 as p —>■ 00 . 

We rescale (13.251) by letting x = $/5, ^ = ^/{SXq) so that 


d ^2 8 1^ 2x ) 

with boundary conditions 

X(l) = Ao/i, x'(l) = pAo/2, 


(3.26) 


(3.27) 


and the far-field condition that x ~ ^o ?/2 as ^ 00 . 

We note that 

X = i>^U - F/7t5)/2 (3.28) 

is an exact solution of (13.261) and, further, satisfies the far-held condition, x 
Aq^/ 2, as ^ 00 . However, this solution is incompatible with the boundary condi¬ 

tion on x^ at ^ = 1, i.e. p = Lq- We conclude that (13.281) is an outer solution with 
an inner, boundary layer-type solution close to ^ = 1 , p = Lq- 
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The inner solution 

The outer solution (13.281) is incompatible with the boundary condition on x'(l); 
we therefore look for a solution that is a small perturbation to the outer solution, 
but whose derivatives can be large enough to match this inner boundary condition. 
We therefore let 


X = i>^U-F/7r~S)/2 + Ax, f = l+X/~S^/^, 


finding that 


and hence that 


dX2 


2{Xq 


\4 

--Ax 

-F/ttS) 


(3.29) 

(3.30) 


X = ^(Ao'C - + aexp 


2y^{Xl- F/irSy/^ 




. (3.31) 


The boundary condition on x'(l) gives that 


a 


= s-ft 


[xl - fm) 


1/2 


23/2 A, 


o 


-(Ao + F/ttS), 


(3.32) 


i.e. a = 0{S 3/2)_ xhe boundary condition on x(l) then gives, at leading order in 


M = Ao/2 — 


27rAoi5 


(3.33) 


The key quantity of interest is actually Z'{Lo), which we find, from i['(Lo) = 
and 

=Lo + Lo ^l^~^} =0 (3.34) 

P=Lo+ (^O - t) 


dp 


to leading order in 

With the remaining boundary condition at p = Lq determined by (j3.34ll and 
continuity of slope, we now return to the wrinkled region to solve the remainder of 
the problem. 


(d) Return to the wrinkled region 

In ^3E1 we showed that the shape within the wrinkled region was given by 

By considering the outer tensile region, we have now shown that 

C'/di/2 = ^ = Xo/2 - (3.36) 

ZTTAQO 

and that 

Z'{Lo) = 0 (3.37) 
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to leading order in Now, since, by the nature of the boundary layer, the deriva¬ 
tive of a function is larger than the function itself in the boundary layer, we must 
also have Z{Lo) = 0 to leading order, and hence that 


JL = ±- 

27r(5 






while (I3.36P gives us 


and hence 


F Xo — 


M = 


2X1 


3(2 -I- Xq) 

We also recall that from the inner tensile region we had 

2 

= « 1.81527^. 


^2j-2 A/ F/{2'k5) 


XifJ. 




(3.38) 


(3.39) 

(3.40) 


(3.41) 


We have therefore determined relationships for the spring stiffness, F/5^ A/ and ^ 
all in terms of Aq. To close the problem we require an additional condition, which 
comes from considering the displacement across the wrinkled region. 


(e) The final condition 


We recall the relationships for stress and strain in an axisymmetric shell are 
given by Isee IVentsel fc Krauthammer . 2001 . for example) 


and 


dur 


^rr — 


Ur 

^ee = — 
r 



[^drj 


(Trr — verse 

Y 

(3.42) 

1 

fdO 

c 

eree — varr 

(3.43) 

2^ ' 

[dej 

1 ^R~- 

Y 


The dimensionless versions of these equations give results for the horizontal dis¬ 
placement. 

In the wrinkled region we use the condition on to see that 


^ _ C _ 1 ^ 

dp p 2 \ dp) 


(3.44) 


while in the tensile (unwrinkled) regions it is simpler to use that from egg, in 
particular: 

^ (3.45) 

P 

While these expressions hold in different regions, they must agree with the corre¬ 
sponding displacement within the wrinkled region, obtained by integrating (13.441) . 
Now, noting that = 4>'(Lq) = 0 and = 4 /(Lq) = C, eqn (I3.45P 

yields: 

Ur{Lf) + LiZiLJ) = UriL+) + LoZ{L+) = -vC. (3.46) 
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It is therefore convenient to express the continuity of displacement in terms of the 
change in displacement across the wrinkled region, i.e. 


LiZ{Lf) - LoZ{Lo) = [Ur]^^ = / C/; dp. 

JLi 


(3.47) 


We then find that 


48) 


fUo 1 / QZ\ 

LiZ{Lj) - LoZ{Lo) = [Ur]if = Clog{Lo/Li) - ^ + Z dp (3. 

and hence, after integrating the last term in the integral on the RHS by parts: 

Cl„g(Lo/L,) = ))_ -(^) (3.49) 

The condition (13.491) is the final condition required to close the system. It is 
implemented exactly in the numerical results described in the main text. However, 
for the purposes of the theoretical analysis we present here, it is enough to consider 
the scaling behaviour of the different terms. We find that the LHS of (13.4911 ^ 
Ji/2 log 5 while the RHS ~ we therefore conclude that at leading order 


I fdZ 


2 V i9p 


dZ 


o ITT - Pitt ‘^p = O- 


dp 


Letting 


we have that 


f = p/8^/\ ujiO = Z/S 


(3.50) 


(3.51) 


a;(e) = -l + 


2TTpS 


i+e/2-e/idp) 


or 


Zoip) 


= -l + 


F 


27rp5i/2^i/2 6p,\s^/2j 


(3.52) 


Performing the integration over f G (0, Ao) and making use of (13.391) - (13.401) we find 
that 


An + 24An - 36 = 0 


(3.53) 


i.e. 


Ao = 6(\/5-2) Ri 1.I90I3, p = 


0 X 3 F 

a; 0.328944, a; 1.990001. (3.54) 
3(2+ A^) ’5 


To calculate the value of A/ we then return to (I3.41|) from which we find that 

Xi ^ 1.16923 (3.55) 


Supplementary Information 















10 


D. Vella et al. 


(/) The wrinkle amplitude 

We make an ansatz for the perturbations of the shell about the tension-field 
state described above: 


Z(p,e) = Zo(p) -I- /(p) cosmO 


(3.56) 


where /(p) = f{p)-R/ip. In the determination of the energy, the form of the wrinkle 
amplitude, f{r), will be of some interest and so we study this here, making use of 
the results for the mean shell displacement Zo{p), (I3.52|) . 

Substituting the ansatz (13.561) into the stress-strain equations (I3.42D - (I3.43|) we 


find that 

dUr _c 1 fdZoV ^ 
dp p 2 dp y ° 


(3.57) 


and also the “slaving condition” 


m2 p _ C ^ Ur 

4 p2 ^ p p 


(3.58) 


Eqn (13.571) may be integrated once, subject to the condition that Ur{Lo) = —vC — 
Zo{Lo)Lo, to determine Ur{p) within the wrinkled region. Ur may then be elimi¬ 
nated from (I3.58|) to give 


m2 /2 


= -z. 


1^1 Lo 
-Clog — 
P P 


P Jp 


Lo 


Zo + k 


= -C\og^ + - I 
P P P Jp 


Lo 




dZo 
dp 
2^ 


dp 


dp 


(3.59) 


(where we have used Zo(Lo) = 0 to eliminate the constant of integration). 

The integration is readily performed, giving 

m2/2 = (36 - 84 /^) J^log^T 9(7-3^)52 f(l-^2)2^ (3.6O) 


(p) Summary 

In summary we have the asymptotic predictions for the wrinkle positions 

Li ft! I.16923J-1/2 (3 51) 

and 

Lo ft! 1.19013(5^/2 (3 52) 

while the predicted force law is 

F ft 1.99J. (3.63) 

The mean shape of the shell (underlying the wrinkles) in the limit of large 
indentations may be written 

= -1-f i (6- A2)^-piA2^2_ i(2 + A2)y3 Ri -1-p 1.146^+0.708^2-0.854^3 

(3.64) 
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where ^ = p/Lo = p/{X5^/'^). Equation p.64l) is given as (9) in the main text. 
Finally, the slaving condition (13.601) in this limit becomes 


~ 9(7 - 


(3.65) 


4. The sub-dominant energy 

To estimate the number of wrinkles observed, and to understand better the emer¬ 
gence of ‘asymptotic isometry’ in this problem, we consider the energy associated 
with the wrinkle pattern. Unlike the energies 74train and Wgas discussed in the main 
text, this energy vanishes in the limit of vanishing th ickness and so is referred to 
as the ‘sub-dominant’ energy ( Davidovitch et al . 201lll . More precisely, in contrast 
to those “dominant” energies (discussed in the main text), the sub-dominant en¬ 
ergies are proportional to the (square of) wrinkle amplitude /. In principle, there 
are three separate contributions to the sub-dominant energy due to wrinkly undu¬ 
lations: elastic strain (stretching, beyond the energy Z4train that is discussed in the 
main text, that is only associated with the strain in the radial profile), bending 
energy and the excess work done in compressing the gas (beyond 77gas discussed in 
the main text). In practice, the latter energy is proportional to the volume of the 
deformation and so, assuming a sinusoidal wrinkle pattern, there is no additional 
energy stored in the gas as a result of wrinkling. Here, we focus on the former two 
sources of sub-dominant energy. 

In this section we use dimensional variables for convenience. 


(a) Sub-dominant stretching energy 

The requirement that the shear-stress and strain vanish to leading order, are = 
Cre = 0, shows that within the wrinkled region the correction to the radial stress, 
air'^ = Urr — C/r IS linear in the wrinkle amplitude f{r) with 


= pRIp 


C 

r 


Yf 


d^o 

dr^ 


cos m9. 


while the radial displacement may be written 

— /Co cos m9. 


Ur = 


(4.1) 


(4.2) 


This result generali zes an earlier result for a natu rally flat sheet adhered to a stiff 
spherical substrate ( Hohlfeld fc Davidovitch . l2015ll to the case of a naturally curved 
shell. 

To calculate the elastic energy of the system, we must integrate the strain energy 
density, cr^rerr/2 over both r and 9. We are therefore interested only in the parts 
that will persist upon integrating wrt 9 (i.e. those that involve even powers of 
cos m9). We have that 


g — 


dUr 

dr 


= e(0) 
rr 


dr 


c 

R 


i [/' COS m6(] ^ 


d^Co 

dr 2 


/ cos m9 


(4.3) 
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The azimuthally averaged energy density is therefore 


(7 f'T' 


= + i(/y£^ + iy/» (1 - - (4.4) 


In the language of the FT approach to wrinkling ( Davidovitch et al llmil) , the 
dominant energy is the first of the three terms in 631 ) with the other two being the 
contributions to the sub-dominant energy, i.e. the energy associated with wrinkles 
themselves, which vanish with the wrinkles’ amplitude /. The first of these sub¬ 
dominant terms corresponds to the radial stretching of wrinkles, the second to the 
curvature of the wrinkled state relative to its natural configuration. 

We use the dimensionless versions of / and Co to write the two components of 
the strain energy as 


(7'r'r' C'l' 


-(sub) _ 


Y 






(Pz 

dp^ 


(4.5) 


We then have that the stretching energy is given by 
rLo 

Us = TT arr^rr 1' dr 

J Li 

/•I r 


y2 


' Li jLo 


,d^Z 




oc 


V 


Y2 


r3 -2 

-d m , 


(4.7) 


where we have used the slaving condition, (I3.65P . We note that this energy is dom¬ 
inated by the second term on the RHS of (14.61) . 

(6) Bending energy 

The bending energy may be written using the slaving condition for very large 
indentation depths, (13.651) . We find: 

p27z rLo pLo 


Uh = 


oc 


Y 

at leading order in i5. 


pro pro 

\B / r{W^Cf dr de ^ ^B r[m^f{r)/r'^] dr 
J 0 J LI J Lj 

dc 

m'‘6logd (4.8) 


2 T-2 f 

2 ^ Jli,Lo 

BpR 2 


(c) Predicted wrinkle number 


The sum of elastic energies, Wtotai = Ug + Ub is a. function of the number of 
wrinkles, m. This sum is minimized by 


4 

m oc 




2g2 


YB log 5 log S 
and is precisely equation (18) of the main text. 


(4.9) 
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5. Comparison of floating Aims to pressurized shells 


Some of the behaviour exhibited in the wrinkling of a pressurized elastic shell is 
similar to what was seen in a recent study of the indentation of a thin elastic film 
floating on a liquid bath and indented at a point ( Vella et al . 201^. In this section 
we discuss the similarities and difference s between the two scenarios. 

Details of the experiment are given bv IVella et al. ( 201£jl . In brief, thin Polystyrene 
films of Young modulus E = 3.4 GPa, Poisson ratio ly = 0.3, radius = 1.14 cm 
and thickness in the range 85 nm < t < 246 nm were floated at the surface of deion¬ 
ized water, with density pi. In the floating state, the surface tension of the water, 
7 / 1 ,, applies a tension at the boundary of the film and so the stress is initially uni¬ 
form and isotropic, CToo = "fiv The freely floating films were then indented at a point 
with wrinkling beginning at a threshold indentation depth Sc ~ Et)^/"^ where 

= [liv! is the capillary length of the bare liquid-gas interface. As Tabled] 
shows, this behaviour is very similar to what is observed in the pressurized shell ex¬ 
ample with the modifications that in this case the foundation stiffness A'found = Pig 
arises from hydrostatic pressure, and the background tension aoo = 7i«, arises from 
the surface tension of the liquid. 

In the floating film case, the outer edge of the wrinkles very quickly reach 
the edge of the sheet, Lq = Rfiim, at which point the system begins to be in¬ 
fluenced by a wrinkly isometry similar to that observed here. However, since the 
outer edge of the wrinkles is fixed, the length scale over which the film returns 
to be flat, 4urv = is independent of indentation depth. A key qual¬ 

itative difference between the Wrinkly Isometry observed in the floating film and 
the pressurized shell is that in the former there is no unwrinkled isometry (except 
a non-axisymmetric “developable cone” involving delamination of the sheet from 
the liquid sub-phase); in the latter, mirror buckling remains as an allowed, unwrin¬ 
kled, isometry (though, as seen here, it has a higher energetic cost than Wrinkly 
Isometry). 

One interesting feature of the shell problem is that the radial displacement across 
the sheet is not, in fact i5^/4urv, as was found in th e floating sheet problem, and 
might be expected intuitively ('see I Vella et al . 20151 for example). The difference 
between the two problems can be understood as follows: In the floating film problem, 
wrinkles extend to the edge of the sheet and there is no tensile region beyond the 
wrinkles. Because of this, the displacement must be determined by integrating the 
strain in the wrinkled region, which gives a purely geometric result (i.e. one that 
does not depend explicitly on the exerted tension), as expected on the basis of 
asymptotic isometry. For the shell case, there is always an outer tensile region. 
As such, the strain relationships immediately give that Ur{Lo) ^ C ^ The 
difference between the two problems could be attributed to the fact that a naturally 
curved object can accommodate deflection without horizontal displacement because 
of its curvature (just think of the mirror buckling solution in which points merely get 
reflected and so do not move horizontally, even though there is a vertical deflection). 
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Property 

Pressurized Shell 

MB WI 

Floating Membrane 
WI 

^ OO 

pRl2 

■yiv 

Foundation stiffness 

-^found 

- Et/R^ 

pig 

Capillary length 
(aoc/Afound)^/"l 

ip 

4 = {iivlpigf^'^ 

Wrinkle onset 
\Sc-^£c{ao./Et)U^] 

illR 


Lo 

= « i.i%5rY/^ 


•^curv 

0 

0 

4 = 

Li 

increasing with 6 5 


Indentation stiffness 
\k = F/5] 

= TvpR « 1.99pR 

{'IlvPigRtA-m)^^'^ 

u{Lo) — u{Li) 


574urv - 574 


Table 1. Comparison between the two types of isometry of a pressurized shell, i.e. Mirror 
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PACS 46.70.De - Beams, plates and shells 
PACS 62.20. Dc - Elasticity, elastic constants 

Abstract - We consider the point indentation of a pressurized, spherical elastic shell. Previously it 
was shown that such shells wrinkle once the indentation reaches a threshold value. Here, we study 
the behaviour of this system beyond the onset of instability. We show that rather than simply 
approaching the classical ‘mirror-buckled’ shape, the wrinkled shell approaches a new, universal 
shape that reflects a nontrivial type of isometry. For a given indentation depth, this “asymptotic 
isometry”, which is only made possible by wrinkling, is reached in the doubly asymptotic limit of 
weak pressure and vanishing shell thickness. 


Introduction. — Just as one pokes an object to see 
how stiff it is, it is common to measure the mechanical 
properties of small scale objects by indentation, often us¬ 
ing an AFM [T]. In biology and medicine many objects 
of interest are relatively thin capsules swollen by osmotic 
pressure, e.g. viruses [5], bacteria |2], yeast cells [1] and the 
polymerosomes used in drug delivery Recent studies 
have focussed on modelling the indentation of pressurized 
elastic capsules [5]|31IS1IB] using numerous approaches in¬ 
cluding finite element methods and analytical techniques. 

The problem of the indentation of a shell with no in¬ 
ternal pressure (pressureless) has been extensively studied 
HHU. In this case, the classical picture is that a spher¬ 
ical shell of thickness t and radius R will accommodate 
indentation by ‘mirror-buckling’: it is as if a portion of 
the shell were sliced off, inverted and then reattached to 
the remainder of the shell (fig. [I^). In reality, the appar¬ 
ent discontinuity in slope is resolved by a boundary layer 
of width (ti?)i/2 [TU]. Since mirror-buckling represents a 
simple isometry of the shell (i.e. there is no strain except 
at isolated points or along curves), the energy of deforma¬ 
tion is localized within this boundary layer and vanishes 
in the limit t —>■ 0. This simple geometrical picture has 
recently been shown to hold even in the presence of an 
internal pressure, p, albeit under the assumption that the 
shell remains very close to an axisymmetric shape [B]. 

However, it is well known that deformed shells often do 
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Fig. 1: (a) Schematic illustration of the mirror buckling of a 
pressurized shell: an indentation depth S is accommodated by 
the (energetically cheap) formation of an inverted cap of width 
{SRyO^ (p) jji reality, an indented pressurized shell wrinkles 
[7]. (c) Setup and notation for the model calculation of point 
indention of a pressurized spherical shell presented here. 


not remain perfectly axisymmetric. Under point indenta¬ 
tion, pressureless shells develop vertices [I2UI6] while pres¬ 
surized shells wrinkle (fig. Hb) [7]. A similar phenomenol¬ 
ogy is observed in other deformation modes including the 
contact between a plate and a shell [HI HZ] and the sec¬ 
ondary buckling of depressurized shells mils]. 

As was shown in previous work [7], the appearance of 
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wrinkles upon the point indentation of a pressurized shell 
is indicative of azimuthal compression. Since the mirror- 
buckled shape remains an isometry of the shell, it is tempt¬ 
ing to assume that the wrinkled shape will be a small per¬ 
turbation of this shape, as has been proposed previously 
[7]. However, here we show that wrinkling changes the 
shape of the shell markedly. Our analysis reveals a new 
wrinkly isometry of the shell that is energetically cheaper 
than mirror buckling. 

Numerical experiments. — Numerical simulations 
were performed using the commercial finite element pack¬ 
age ABAQUS (SIMULIA, Providence, RI) for a spher¬ 
ical shell of radius i? = 1 m, thickness in the range 
0.3 mm < t < 2 mm. Young’s modulus E = 70 GPa, and 
Poisson ratio ly = 0.3. The shell is subject to an internal 
pressure p — 0(10® Pa). (Four-node thin shell elements 
with reduced integration were used in all calculations and 
a mesh sensitivity study was carried out to ensure that the 
results are minimally sensitive to the element size. No ini¬ 
tial geometrical imperfection was used and the ABAQUS 
automatic stabilization scheme, with a dissipated energy 
fraction of 2 x 10“^ and an adaptive stabilization ratio of 
0.05, was employed to capture instabilities in the shell.) 

Simulations show that for sufficiently small indentation 
depths (5, the shell remains perfectly axisymmetric. In this 
limit the (thickness integrated) in-plane stress [11] is close 
to the isotropic tension prior to indentation: CToo = pR/2. 
Deformation of the shell is resisted by a normal force that 
is proportional to the vertical displacement — effectively 
a Winkler foundation of stiffness [5|. This founda¬ 
tion stiffness depends on the stretching stiffness of the 
shell, Et, and so on dimensional grounds we might expect 
that Kw ^ Et /an analysis of the governing equations 
confirms this expectation BE- By analogy with a liq¬ 
uid’s meniscus, the combination of an intrinsic tension, 
O’oo ^ pR, and a linear foundation stiffness, ^ Et/R?, 
gives rise to a characteristic horizontal length over which 
vertical deflections decay 


P — 

Up — 


J V / 


( 1 ) 



Fig. 2: The evolution of the inner and outer wrinkle posi¬ 
tions, Li and Lo, with increasing indentation depth S. The 
results of ABAQUS simulations are shown for different values 
of the dimensionless pressure r and thickness ratio t/R: t = 
26.2, t/A = 3 X 10“^ (red triangles), r = 52.5, t/A = 3 x 10“^ 
(green squares), r = 177,1/R = 2 x 10“® (blue circles) and 
T = 378, t/R = 1 X 10~^ (magenta diamonds). The predictions 
of the FT membrane theory (valid as r ^ oo) are shown (solid 
black curve) together with the asymptotic results (HOI) and H3II 
(dashed lines). For comparison, the results of the axisymmetric 
membrane-shell theory (3 are also shown (dotted curve). 


the shell, it is also similar to the “tensional bendability” 
introduced in recent studies of uniaxially stretched sheets 
m- We neglect the effects of bending since all simula¬ 
tions reported here have r > 25. In the language of [21], 
we consider the limit of large bendability. 

The picture outlined above holds while the stress within 
the shell remains close to its pre-indentation level. How¬ 
ever, indentation also introduces a geometrical strain e ^ 
{6/£p)^, and a corresponding contribution to the stress 
~ Et5^/ip. This additional stress leads to an inhomogene¬ 
ity in the stress within the shell: the radial stress arr/r) 
decreases monotonically to the far-held value, CToo = pR/"2, 
while the hoop stress agg decreases more rapidly and over¬ 
shoots the far-held value (since material circles are pulled 
inwards and become relatively compressed) [7j|22]. The 
relative size of the geometrical stress ^ EtS'^/i'^ and the 
pre-stress ^ pR is a function only of the ratio: 


The length ip is an effective capillary length. 

Another restoring force is due to the bending rigidity, 
B = Et^/[12(1 — z/^)]. However, previous work [5)17] has 
shown that the effect of bending may be neglected when 
the internal pressure, p, is sufficiently large. To see this, 
we note that the foundation stihness resulting from bend¬ 
ing effects ^ B/L^, which is signihcantly smaller than the 
foundation stiffness ^ Et/R/^ (discussed above) when 


pR/^ 

jEtB)^ 


» 1 . 


( 2 ) 


The parameter r may be thought of as the dimensionless 
pressure within the shell. Since r is the ratio between 
the tension-induced and bending-induced stiffnesses of 


S = SR/il = EtS/pR^. (3) 

Here, this dimensionless indentation depth plays the role 
of the “confinement” parameter introduced in previous 
studies [221 - 124] . Above a critical value, d > dc ~ 2.52, 
the hoop stress becomes compressive, agg < 0, beginning 
a distance Tc ~ 0.58£p from the indenter and confined 
to an annulus Lj < r < Lq as 6 increases [7]. Since 
we assume that r ^ 1, (the shell has negligible resis¬ 
tance to compression), any compressive stress immediately 
leads to wrinkling. Assuming that the shell remains close 
to its axisymmetric state, detailed calculations [7] show 
that as S increases, wrinkles are confined to an annulus 
Li < r < Lo with both Lj and Lq increasing (see fig. [2]) . 

Simulations confirm the behaviour at the onset of wrin- 
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kling [7], in particular the prediction Sc ~ 2.52. However, 
as fig. [5] shows, systematic discrepancies are observed in 
the wrinkle position as indentation increases beyond Sc- 
In this letter, we show that this discrepancy signals an 
important change in the mechanical response. 

Tension field theory description of wrinkling. — 

The key ingredient missing from previous works IS [7] is 
that the hoop stress, age, cannot become arbitrarily large 
and negative. Instead, the appearance of wrinkles allows 
the shell to relax the compressive stress so that agg Ri 0, a 
result that is referred to variously as Tension field theory 
[25ll^ . Relaxed energy [27ll^ or Far from threshold (FT) 
wrinkling [20112^. 

An immediate consequence of the collapse of hoop 
stress, agg ~ 0, is that the radial stress arr ~ C/r in 
the wrinkled region {Lj < r < Lq), for some constant C. 
We denote the normal displacement of the wrinkled shell 
from its natural state by ^(r, 9) = ^(r) + /(r) cos m6: ( is 
the mean displacement upon which m fine wrinkles of am¬ 
plitude /(r) are placed. (Here r is the horizontal distance 
from the axis of indentation and 6 the azimuthal coordi¬ 
nate, see fig.Ht.) In the FT approach, ({r) is determined 
at leading order with m and / determined subsequently 
[20l[23]. In the wrinkled zone vertical force balance reads 

m 

which can immediately be integrated to give 


c 


^2 3 p 

2R 6C 2ttC 


(5) 


for some constant of integration A and the force F. 

To determine A, F, and (7, together with the inner and 
outer limits of the wrinkled region, Lj and Lq, the above 
solution must be patched together with the solution of 
the vertical force balance in the tensile (i.e. unwrinkled) 
regions 0 < r < Lj and r > Lq, which reads [71I5D] 


lA. (L 

r dr \R 



= P- 


F 6{r) 
2TT r 


( 6 ) 


Here i/jir) is the Airy potential, such that a^-r = 4’/f 
agg = which satisfies the compatibility condition [7] 


1 A 

Y dr 


Id 

(r^/>) 

r dr 


1 

i? dr ^ \ dr / 


(7) 


Each of the equations in the tensile regions, m and m, 
is second order; there are therefore 4 degrees of freedom 
associated with each tensile region. In the wrinkled region 
there are 3 degrees of freedom {A, C and F). Finally, with 
the radii of the wrinkled zone L/, Lq, we have 13 degrees of 
freedom in total, which must match the number of bound¬ 
ary conditions (BCs). At the indenter, r = 0, we impose 
the vertical displacement, C(0) = —d, while the radial dis¬ 
placement vanishes, rir(O) = 0. Far from the indenter, we 


require that C(r) —>■ 0, while ip ~ pRr/2, the uniformly 
stressed solution. In addition, we also require continuity 
of C, 4', arr, o'ee and Ur, at r = Lj and r = Lq- The 
first four of these continuity conditions at each boundary, 
together with the conditions at 0 and oo, give a total of 
12 BCs. Imposing continuity of Ur at Lj and Lq requires 
the integration of a further equation for durjdr across the 
wrinkled region, introducing a new unknown and two BCs. 
After some algebra a single condition emerges 


C, Lo 

y‘“®L7 



r dC 


( 8 ) 


Here, the LHS and RHS represent, respectively, the 
stretching of lines of longitude, and how the boundary 
displacement differs from that expected were the shell to 
deform inextensibly. 

Equations can be solved numerically using 

the MATLAB multipoint boundary value problem solver 
bvp4c. For a given <5 > (5c, the solution gives the inner 
and outer limits of the wrinkled region (L/ and Lq re¬ 
spectively), together with the force F. Figure [5] shows 
the evolution of L/ and Lq with (5, together with the re¬ 
sults of our ABAQUS simulations. Also shown for com¬ 
parison are the extent of the compressive hoop stress for 
a perfectly axisymmetric membrane-shell [7]. The agree¬ 
ment between simulations and the results of this collapsed- 
compressive stress model, which takes into account the 
non-perturbative effect of wrinkling on the shell shape, is 
significantly better than that based on the axisymmetric 
membrane-shell theory. In particular, we observe that L/ 
decreases with increasing b — qualitatively different from 
the predictions of a standard post-buckling analysis [7]. 

Figure |3] shows the force-displacement relationship: 
ABAQUS simulations and the numerical solution of ([S])- 
(O produce results that are in good agreement with each 
other. Furthermore, in the limit of large indentation 
depths, the force-displacement relationship becomes lin¬ 
ear, corresponding to a constant stiffness k = F/pRS (in¬ 
set of fig. ED- At the scaling level, this result is in qual¬ 
itative agreement with prior analysis of the axisymmet¬ 
ric membrane-shell equations and with the simple picture 
of deformation being accommodated via mirror-buckling 

[5] : the work done ~ FS is stored purely in the gas com¬ 
pressed at constant pressure pAV ^ p[{6Ry^‘^]‘^5, so that 
F ~ pRb. However, for the axisymmetric, mirror-buckled 
shape, a detailed calculation predicts that k = F/pR6 = tt 

[6] while fig. |3] indicates that k ~ 2. To gain quantitative 

understanding of the differences between the axisymmet¬ 
ric membrane theory and the compressed-stress approach 
described so far, we now consider the limit of very large 
indentation depths, (5 1. 


Analytical results for large indentation. — As 

already remarked, L/ —>■ 0 as (5 increases. This can be 
shown by a more formal calculation [30] but for now we 
make this an ansatz that we verify a posteriori. With 
this assumption, it is clear that the constant A in ([5]) is 
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Fig. 3: Numerical results for the force-displacement relation¬ 
ship of an indented pressurized shell. The results of ABAQUS 
simulations are shown for a shell thickness t/R = 3 x 10““*. 
Different pressurisations are used: p = 0.1 MPa, r ~ 52.5 (red 
circles) and p = 0.05 MPa, r « 26.2 (blue triangles). The nu¬ 
merical solutions of are shown (solid black curve) with 

the asymptotic relationship (1121) (dashed black line). The cor¬ 
responding results for a perfectly axisymmetric shell are also 
shown (dash-dotted curve) with the ultimate mirror-buckling 
result F = nSpR (dotted line) [5]. The black square shows the 
onset of wrinkling, at <5 « 2.52. 


A —S = C(0)> the imposed indentation depth. A de¬ 
tailed calculation m shows that outside the wrinkled re¬ 
gion, r > Lo, the vertical deflection of the shell decays 
exponentially. For the wrinkled zone to patch onto this 
we must have C(^o) = C'[Lo) = 0. The deflection of the 
wrinkled portion of the shell may then be rewritten 

J=-l+(f-^)f + Y«^-(2 + A’)«V4 (9) 

where ^ = rjLo and A = Lo/{6Ry/‘^. Substituting ([9|) 
into the constraint (0), and neglecting the logarithmic 
term, we find A = [6(\/5 — 2)]^/^ Ri 1.19 so that for 5 1: 

Lo Ri 1.19013((5i?)i/2 (10) 

Comparing coefficients between Eqs. (0) and (0) we find: 

C = 3(1 0.329pS^/^R^/^ , (11) 

To determine how L/ vanishes with increasing inden¬ 
tation depth we use a scaling argument. Recall that the 
radial stress in the wrinkled region is given by arr = Cjr. 
The stress in the inner tensile region (r < Lj) may be es¬ 
timated by using the geometrically induced strain in the 
deformed region, r < Lq, giving arr Et{5/LoY 
similar estimate for the stress in a tensile inner core has 
been used in a related problem m)- Equating these two 


stresses at the inner edge of the wrinkles, we find that 
Lj ^ pi?®/^(5“^/^/(Et). A detailed calculation [SD] yields: 


Li 


1.16923 



(13) 


The predictions (fTUl) . (fT^ and (fT^ are in excellent agree¬ 
ment with the solution of the full tension-field problem 
and with ABAQUS simulations (figs 0 and 0). 

We note that the force law in m has the same scal¬ 
ing behaviour as the perfectly axisymmetric case [B], but 
with prefactor 1.99 rather than tt. In particular, this scal¬ 
ing is independent of the elastic properties of the shell {E 
and t) depending only on the internal pressure and shell 
radius. In the perfectly axisymmetric shell, this indepen¬ 
dence of material properties was a natural consequence 
of mirror-buckling being an isometry of the shell 0. We 
shall see that an asymptotic isometry of the wrinkled shell 
(i.e. with the compression collapsing faster than the resid¬ 
ual tension), underlies the form of the force law in (00) . 

Finally, using the various pre-factors determined by the 
preceding analysis, we rewite Eq. 


c 

(5 


-I -b 0.963 


r 

(5i?)l/2 


1 

2AR“'^'^°^(5i?)3/2’ 


(14) 


for 1.17^p/((Ii?)3/2 < r < I.19(<5i?)3/2. When written in 
this way the profile in (fHll appears universal. Fig.H]shows 
that cross-sections of simulated shells tend to this univer¬ 
sal shape as the confinement 5 increases. In particular, 
we see that (TH)! gives a much better account of numeri¬ 
cal simulations than classic mirror-buckling, which, in this 
notation would read (^/S ~ — I -b / [5R). A puzzling ob¬ 
servation is the asymmetry of wrinkle crests and troughs: 
the profile (11^1) captures the former better than the latter. 

The energy stored in the indented shell. — Let us 

evaluate now the elastic energy stored in the shell at large 
indentation depth, when L//£p —^ 0. This energy consists 
of two principal parts: the stretching energy ^strain stored 
in the tensile component arr = C/r of the compression- 
free stress field, and the bending energy 74>end of the wrin¬ 
kles. The evaluation of ^strain is straightforward: 


^strain = j Tarr^rr dr -^l51og(5 , (15) 


where we used Eq. m- A more careful calculation of 
Z^strain, which accounts also for the tensile core zone r < 
Lj, confirms this scaling relation. 

A calculation of Uhend, similarly to previous works m 
[23], proceeds in two steps. First, we express the energy as 
a function of m only, by employing a “slaving” condition 
between the number, m, and amplitude, /(r), of wrinkles: 


2 f 2 
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R dr 
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Fig. 4: Approach to the wrinkly isometry of a pressurized shell. 
Cross-sections through an indented pressurized shell, obtained 
from ABAQUS simulations, are shown for different indenta¬ 
tion depths with r « 52.5: 5 = 6.1 (red circles) and 5 = 15 
(green squares); to reach larger indentations numerically we 
must consider lower pressure, r = 26.2 and S ~ 90.3 (blue tri¬ 
angles). Cross-sections through a trough in the wrinkle pattern 
are shown on the left (r < 0); Cross-sections through a crest 
are shown on the right (r > 0). The predictions of the FT the¬ 
ory with the corresponding (finite) value of 5 are shown by the 
coloured dash-dotted curves. As 5 increases the deformation 
tends to the wrinkly isometry (1141) (solid black curve), rather 
than the classic mirror-buckled isometry (dotted curve). 


Here the RHS may be calculated from the results of the 
preceding analysis, Eqs. (I5I10I11I12I) . Physically, (fTBl) 
stems from the fact that, for the compressive stress to 
collapse, the azimuthal undulations must accommodate 
the excess length of deformed latitudes. (Notice also that 
Eq. ([5|) can be obtained from the slaving condition (fTCl) 
since the wrinkle amplitude vanishes at r = L/, Lq-) 

In the second step, one determines the energetically fa¬ 
vorable value of m, for which the bending energy balances 
the contribution of wrinkles to the stretching of longitudes 
(beyond ^strain)- Intuitively, bending favors small curva¬ 
ture (i.e. small m) while stretching resists deviations from 
C{r), thus favoring small /(r) or, by condition (ITCl) . large 


m. The detailed calculation [5^ yields: 


\\og5j 

(17) 

and gives the bending energy 


Wbend ^ pi?h^(log5)^/^/T. 

(18) 


As shown above, our ABAQUS simulations are in excel¬ 
lent agreement with the macroscale features of the wrinkle 
pattern predicted theoretically (L/, Lq, F and the shell 
profile). However, while close to threshold (i.e. 5 = 0(1)), 
previous ABAQUS simulations [7] confirm the r depen¬ 
dence predicted by (EZl), the predicted evolution of m with 
increasing 5 has not been observed. This shortcoming is 
consistent with recent work, which showed significant hur¬ 
dles in numerically predicting the wrinkle number [31] . In 


that study, which addressed a planar sheet under differ¬ 
ential tensile loads, it was shown that under-damped dy¬ 
namic relaxation allows “hopping” between meta-stable 
states, leading ultimately to the ground state (i.e. the 
wrinkle number of the global energy minimizer). In our 
simulations, the lack of variation of wrinkle number per¬ 
sists despite our use of an over-damped relaxation tech¬ 
nique, suggesting either a sensitivity to numerical tech¬ 
nique or a higher energetic barrier between meta-stable 
states than in other problems [31]. Additionally, a more 
elaborate calculation of the energy may result in a spatial 
profile of the wrinkle number, m —>■ m(r) pSl 15^155] . 

Wrinkly isometry. — A key feature of (11^1) in the 
limit of very large indentations, 5 1, is that the wrinkle 

pattern covers almost the entire region from the indenter 
to the outer ridge, which is positioned at r Ri 1.19(5i?)^/^; 
outside this ridge, the shell is effectively undeformed |5D] . 
As such, the work done in compressing the gas is: 

Ugas = pAU = 27rp f rC dr Ki 0.995pRS^. (19) 

Jl, 

Interestingly, the derivative dWgas/d(5 is identical to the 
total force F calculated above (Eq. [12]). The origin of 
this remarkable feature is elucidated by comparing Ug^s to 
the energies of stretching and bending stored in the poked 
shell. We have that Ustram/Ugas ^ (log^)/(5 <C 1, since 
(I14|l only holds for 5 1. At the same time, for the FT 

expansion [20] to be valid we must have Uhend ^ Ugas, 
i.e. T (logh)^/^. The two components of the elastic 
energy are therefore negligible in comparison to the work 
done in compressing the gas at constant pressure; this is 
the signature of the “asymptotic isometry” that is realized 
by the wrinkled shell. 

Rather than a specific shape, such as mirror-buckling, 
which exists also in the absence of pressure, and to which 
the undeformed shell can be mapped isometrically almost 
everywhere (except at the circular ridge), a wrinkled shell 
develops finer and finer structure in a macroscopic zone 
(L/ < r < Lo)- However, this isometric shape can only 
be approached when both Z4>end,ZYstrain ^ ^gas, he. when 

1 < 5 exp(r^). (20) 

We shall show shortly that this occurs when both the shell 
thickness and the pressure become asymptotically small. 
Such a doubly asymptotic approach to isometry has been 
noted recently in set-ups that involve strong deformation 
of thin sheets; in these systems the work done by the ex¬ 
ternal load is generally stored elsewhere than the deformed 
sheet itself pTlI^IM] . 

Conclusion. — Why is the wrinkly isometry selected 
by the poked shell, and not the simpler, mirror-buckled 
isometry? In the parameter regime (1201) the deformation 
energy is well-approximated by Wgas alone: the energies of 
bending and stretching (which are localized in the circu¬ 
lar ridge of the mirror-buckled shape but spread around in 
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the wrinkly isometry) are insignificant. The only energy 
that one should compare then is the work of compressing 
the gas. In both types of isometric response this energy 
oc pRS^ but for wrinkly isometry the prefactor ft! 0.995 
while for mirror-buckling it is tt/2: since wrinkly isome¬ 
try compresses less gas than mirror-buckling, it is, in this 
limit, energetically favourable to mirror-buckling. 

To understand the physical nature of the asymptotic pa¬ 
rameter regime (EOI), we note that 0 may be re-written 
1 S /{ tT ) <C exp(r^). For a fixed indentation depth 
(5, the first inequality may be achieved if, tr —^ 0; how¬ 
ever, the second inequality then holds only if t —>■ 00. If 
the Young’s modulus and radius of the shell are consid¬ 
ered fixed (as in our simulations) this may be achieved 
with vanishing thickness and pressure, t,p ^ 0, but with 
the additional requirements that p/t — >■ 0 and pjt ^ —^ 00 
(e.g. p ^ for any 1 < a < 2). The wrinkly isometry 
therefore exists if the pressure is sufficiently weak, but not 
too weak. 

If instead, p/t^ —>■ 0 as t —>■ 0, the pressure is “ultra- 
weak” . In this case, the energy of deformation is not well- 
approximated by Ugas, and the bending energy of defor¬ 
mation must be considered. One might guess that in such 
a case the deformed shell will approach isometry simply 
by selecting the mirror buckled shape. A host of experi¬ 
mental and computational studies mill], as well as one’s 
own experience of poking a ping-pong ball, indicate that 
such a guess is naive: instead, a polygonal shape is formed, 
with a number of vertices that increases slowly upon in¬ 
creasing indentation depth (relative to thickness). Such 
a shape is strictly different from both the axisymmetric 
mirror-buckled shape and the wrinkly isometry described 
here (for a weak, but not ultra-weak, pressure). This sug¬ 
gests the fascinating possibility that a new type of isome¬ 
try characterizes the poked shell for small thicknesses and 
ultra-weak pressure. We hope that future work will clarify 
the nature of this third type of isometry of shells. 
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